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Abstract 

Recently Wang et al. carried out a laboratory experiment, where a Brownian particle was 
dragged through a fluid by a harmonic force with constant velocity of its center. This experiment 
confirmed a theoretically predicted work related integrated (I) Transient Fluctuation Theorem 
(ITFT) , which gives an expression for the ratio for the probability to find positive or negative values 
for the fiuctuations of the total work done on the system in a given time in a transient state. The 
corresponding integrated stationary state fluctuation theorem (ISSFT) was not observed. Using 
an overdamped Langevin equation and an arbitrary motion for the center of the harmonic force, 
all quantities of interest for these theorems and the corresponding non-integrated ones (TFT and 
SSFT, resp.) are theoretically explicitly obtained in this paper. While the (I)TFT is satisfied for 
all times, the (I)SSFT only holds asymptotically in time. Suggestions for further experiments with 
arbitrary velocity of the harmonic force and in which also the ISSFT could be observed, are given. 
In addition, a non-trivial long-time relation between the ITFT and the ISSFT was discovered, 
which could be observed experimentally, especially in the case of a resonant circular motion of the 
center of the harmonic force. 



1 



I. INTRODUCTION 



Fluctuations of physical properties of statistical mechanical systems were first considered, 
in the modern context of dynamical Hamiltonian or dissipative systems theory, by Evans, 
Cohen and Morriss[0]. It concerned here the statistics of phase space contraction or entropy 
production fluctuations over a certain time interval. In particular, the probabilities for equal 
positive or negative entropy production fluctuations of a certain magnitude were considered. 
Two different physical situations have been treated. First, in Ref. for a non-equilibrium 
stationary state, possibly far from equilibrium, the fluctuations of the dissipative (viscous) 
part of the pressure tensor of a fluid were studied. Next, Evans and Searles0 studied the 
fluctuations of entropy production in an ensemble of phase space trajectories emanating from 
an initial equilibrium state in the course of time. While the first case concerned a study 
of stationary state fluctuations in trajectory segments of a given duration along a single 
trajectory in a non-equilibrium stationary state and will be called the Stationary State 
Fluctuation Theorem (SSFT), the second case involved a study of an ensemble of many 
transient phase space trajectories each over a time r, all emanating from an equilibrium 
ensemble at time t = 0, which will be called a Transient Fluctuation Theorem (TFT). 

Mathematical proofs have been given of both theorems 0, ^, ^ ^ and many computer 
simulations have confirmed both theorems(e.g. [jl], 0, |^). While the original proofs of both 
FT's were based on the deterministic dynamics of many particles, later proofs for systems 
with stochastic dynamics were given by Kurchan]^ and Lebowitz and Spohn0. Only one 
laboratory experiment had been carried out for the SSFTp and none for the TFT, until 
recently by Wang et al.^. 

All deterministic theories were concerned with systems in phase space consisting of many 
particles. The experiment of Wang et al. was carried out for a single Brownian particle 
which was dragged by means of a uniformly moving harmonic potential generated by a laser 
through a many particle molecular solvent. This system differs in an important aspect from 
the many particle systems in the phase space considerations. This is due to the fact that 
the (mesoscopic) Brownian particle is much heavier than the surrounding fluid particles, 
which makes it tractable in a different, much simpler though approximate, way from the 
full dynamical systems treatment in phase space mentioned above. In fact, the treatment 
generally applied to such systems is via a Langevin equation for the stochastic motion of 
the Brownian particle in a medium in real space, which is characterized only by its friction 
with the particle and its temperature. As a consequence the very complicated many particle 
problem can be treated by a single particle Langevin equation, if it is near equilibrium and 
on the level of Irreversible Thermodynamics [|10|] . 

However, for the investigation of the Fluctuation Theorems, an additional difficulty is 
that the experiment considered here, uses a time dependent force on the particle, since 
the Langevin equation contains a laser-induced harmonic force on the particle, where the 
position of the minimum of the harmonic potential changes in time. As a consequence, the 
treatments in Refs. 0] and do not directly apply to this experiment. 

Furthermore, the phase space treatments of dynamical systems have always been such 
that (at least if the total energy of the system is kept constant) the total phase space 
contraction can be directly related to the total entropy production of the system. This 
has led to the TFT0, Q and the SSFT0, ^ ^ for the entropy production. To be sure, 
this connection between phase space contraction and entropy production can only be made 
if the total work done on the system is purely dissipative. However, in the Wang et al. 
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experiment, this is not so. 

To see this, it is useful to consider the total (tot) work W^"^ done on the system during 
a time r: 

iy*°*= rrftv;-F(x„x,*), (1) 

^0 

where 

F(xi,x:) = -Mxi-x:). (2) 

is the harmonic force exerted on the particle, with x^ the position of the particle and xj' 
the position of the minimum of the harmonic potential, and k the force constant of this 
potential. Furthermore, in Eq. (|l|) = xj'. At t < 0, in the Wang experiment, the center 
of the harmonic potential is at rest at Xq = 0. At t = 0, the harmonic potential is set in 
motion relative to the fluid with a constant velocity v*, so that x* = v*t for t > 0||11||. 

The crucial question is now what is the dissipative part of which is responsible for 
the heat or entropy produced in the system as a result of the friction of the particle with 
the surrounding fluid. The dissipative part should not include any purely mechanical work. 
To see how W^°^ is related to the dissipated work over a time r, we rewrite the total work 
done, W^°^ in Eq. (|1|), as follows 

= fdt^ji-n^t.^i) 

Jo 

= - dt{wt- V*) ■ F(xi, x,*) + / dtwf F(xi, x^) 

= k i dt{kt- k;) ■ (xt - X*) + / dtvf F(xt, X*) 
Jo Jo 

= AU + W^\ (3) 
defining AU = | [|Ax^p - |Axop] with Ax^ = x^ - x*, and 

iyf^= ^dt^rt■F{^t,^;)■ (4) 



Here, W^^ is the work done on the Brownian (Br) particle by the harmonic force. As (at 
least ideally) the Brownian particle has no internal energy, all this work is converted into 
heat, which is the source of the entropy production. Hence, W^^ is the dissipated work. On 
the other hand, the term AU in Eq. represents the purely mechanical ( "center of mass" ) 
work done on the particle in the external harmonic potential. 

Therefore, Wang et a/.'s entropy production during time r, denoted by S,- in Ref. is 
really the total dimensionless work done on the system, and we will denoted it by 



Wr = I3WI°' = (3 / rftv* ■ F(x„x*), (5) 
Jo 

where /3 = l/lksT), with ks Boltzmann's constant and T the temperature of the surround- 
ing fluid. By following the position of many, independent Brownian particles and using 
Eqs. (|^) and (|^), Wang et al. measured this dimensionless work Wr — or what they called 
the entropy production H-r — over time intervals r and constructed from that the probability 
distribution function P{Wr), which they found satisfies 

^(^-) = (6) 
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So they established experimentally the validity of a TFT for the total work done on the 
system {f3W^°^), rather than for the entropy production of the system. Strictly speaking, 
Wang et al. measured an integrated variant (an ITFT) of the TFT in Eq. (^), explained in 
Sec. [11 D| . A direct transformation of the TFT for Wr to a TFT for the dimensionless 
entropy production, which would be jSW^"^ , is not obvious, since At/ in Eq. (^) fluctuates 
also. 

We remark that Mazonka and Jarzynski[|T2| studied the same system as used in the 
experiment of Wang et al. theoretically — before the experiment — and derived the TFT 
and the SSFT for the total work done on the system, but not for the entropy production |]T3 



Unaware of Mazonka and Jarzynski's work, but in view of the experiment of Wang et 
al., we studied this experiment independently [jl4|. For, Wang et a/.'s experiment is clearly 
important for practical purposes, since it involves a general property of the work done on 
a system. We discuss the observability of the work related TFT as well as SSFT for an 
arbitrary rather than a uniform motion of the harmonic potential. So, for the purpose of 
treating the experiment and its generalizations, in this paper we too will treat the TFT and 
SSFT for the dimensionless work, with a focus on the feasibility to do a convincing SSFT 
experiment. To the best of our knowledge, no Fluctuation Theorem for entropy production, 
either an (I)TFT or an (I)SSFT, has been derived for a Wang-type system, neither from a 
phase space perspective nor in real space (via a Langevin equation). 

The outline of the paper is as follows. In Sec. |I| we present our Langevin model and 
we develop the general theory for the verification and the experimental observability of the 
work related fluctuation theorems, and discuss a new relation between the fluctuations in the 
transient and in the stationary state. In Sec. we specialize the general theory to the case 
of a linear and a circular motion of the minimum of the harmonic potential, investigating in 
detail the observability of the ITFT and the ISSFT. In Sec. |V| we end with a discussion. 



II. THEORY 



A. Definition of the Model 

Like in the experiment of Wang, the model we consider has a spherical Brownian particle 
in three dimensions with a radius R and mass m in a fluid with viscosity r] and temperature T 
and the Brownian particle is subject to an external harmonic potential with a time dependent 
position of its minimum. For t < the minimum of the harmonic potential is at the 
origin, xj' = 0, whereas for t > it moves with a velocity v^, which can be, in principle, 
an arbitrary function of time. The equations of motion for the particle are then of the 
Langevin-type: 

±t = Vi (7a) 

m\t = -awt - k{x.t - x.^) + (7b) 

where and are the position and velocity of the Brownian particle, respectively. In this 
equation, the Brownian particle feels three forces. The first force is the drag force — av^, 
with according to Stokes' law 

a = 67ir]R. (8) 
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The second force is due to the harmonic potential [see Eq. (^]. The third and last force is 
a random force (t, which is taken to be Gaussian and delta-correlated in time: 



iCt) = 0; iCtCs) = 2kBTa5{t - s) 



(9) 



The strength of the random force in Eq. is such that the equilibrium distribution function 
Peg for X and v 



eq 



27T 



(10) 



is stationary under the equations of motion Eqs.(^b)[|T^. 

The system will only be considered in the strongly overdamped case 

mk < (11) 

Effectively therefore, the mass can be seen as a small parameter and will be set equal to 



zero 



15|. From Eqs. (|7a]) and (7B), we find then a simplified Langevin equation for the 
position of the particle only 



7:-i(x,-x*)+a-iC 



with a relaxation time 



a 

= k- 



(12) 



(13) 



(14) 



It is convenient to separate the average motion of the Brownian particle (which results 
from the deterministic forces alone), from the stochastic motion. The average motion is 
given by the solution yj' of the deterministic part of the Langevin equation (^2|), i.e., by 



When we only use Xf, the equilibrium distribution in Eq. (|T^) reduces to 

Pe,(x) = [ dwf,^{^,^) = (fc/3/27r)3/2e-/^tl-l^ 



(15) 



with initial condition Yq = 0. We can then look at the deviations from this average motion 
by introducing the transformation 



X< = xt - y;. 

This turns the Langevin equation (IT^) into the simple form 

±t = -r-'Xt + a~\t. 
yl follows from the general solution of Eq. (15): 



Jo 



(16) 



(17) 



(18) 
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so that with Yq = 0, and a partial integration, one obtains 



*dt'e-(*-*')/"'-v;,. (19) 







The transformation (|T6D with Eq. (|l^) can be interpreted as going to a co-moving frame, 
but it is not co-moving with the minimum of the harmonic potential, but with which is 
what the motion of a particle starting at Xq = would be if there would be no noise term 
in the Langevin equation (0). 

Eq. ([T7|) shows that in the co-moving frame, one has the standard Ornstein-Uhlenbeck 
process [p!5| , [16|. Its solutions are well-known. The Green's function of the Ornstein- 
Uhlenbeck process, which gives the probability for the particle to be at Xi at time ti, 
given that it was at Xq at time to, is Gaussian in both Xq and Xi. Its stationary solution is 
of the form peg(X), with peg given in Eq. ([T^). Initially the particle is distributed according 



to Eq. (0), but because Xq = xq (yq = 0), one sees that the initial distribution is already 
the stationary one, and in this special, co-moving coordinate frame, the distribution of the 
Brownian particle has an equilibrium distribution for all time: 

P(X,t) = (/3A:/27r)3/V^II^I'. (20) 

We end this section by writing Wr in Eq. (|^) in terms of X^, 

= -k(3 fdt [v* • X, + v,* ■ (y; - x,*)] . (21) 

^0 



B. Transient Fluctuation Theorem for the Total Work 



In Eq. (|2lD , W-r is a linear function of X^. Combined with the Gaussian nature both of 
the Green's function of the Ornstein-Uhlenbeck process [Eq. ([T7|) ] and of the initial distri- 
bution[Eq. (|20|) ], this means that the distribution Pt of Wr is Gaussian 



[w-V-Mr(T)]^ 



p 2Vt(t) 

MWr) = , (22) 

^y2■7TVT{T) 

where the subscript T denotes that the transient case is considered. The mean Mt of Wr 
is, from Eq. 

Mrir) = -kp fdt v; ■ (y; - x*), (23) 







since (X^) = [Eq. (^Ol)]. Using the expression for y^ in Eq. ([19D , this can be also written 



as 



Mt{t) = k(3 dt', dt[ e-(*^-*i^/-v4 ■ V,*, (24) 
The variance Vr of Wr is only affected by the first term in Eq. (|2T|) , so that 

Vt{t) = {{Wr-{Wr)r) 

= 2^/3' fdt'^ I 'dt[ v4 ■ (X,, X,, ) ■ v,*, , (25) 
Jo Jo 
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where we used the symmetry of the time-correlation function (X^^Xj'^) under interchange 
of t'l and tg. To calculate this function, notice that Xj has a stationary distribution so it 
can be written as (X^^.^'^Xo). Using the formal solution of the Langevin equation in the 
co-moving frame [Eq. (|I7|)] for t > 0, 

Xt = e-'/^'-Xo + a-' fdt' e-^'-''^'^Xt'. (26) 

one obtains with (Cf) = 0, (Xq) = 0, {CA) = and (XoXq) = [kf3]-% 

(XiXo) = [(3k]-'e-'/^^l. (27) 
The variance in Eq. ([25| ) then becomes 



Vt{t) = 2k(3 dt', dt[ e-(*^-*i)/-v4 • V*, . (28) 

Comparing with the mean in Eq. (p4D , we see 

Vt{t) = 2MT{r). (29) 

This relation leads straightforwardly to the TFT: Given the distribution function of Wr in 
Eq. (p2D, one easily shows that 

^ ' e W , (30) 



Pri-Wr) 
which, by Eq. (PU]), becomes 

P^Wr 



(31) 



Pri-Wr) 

which is identical to the TFT in Eq. (P). 



C. Stationary State Fluctuation Theorem for the Total Work 

To move on to the SSFT, it is necessary to clarify what the stationary state means, 
since in the X coordinate system, the distribution is stationary, which would suggest that 
the TFT is also the SSFT. This is not the case. If one defines a stationary state as that 
state in which (on average) the physical (macroscopic) parameters do not change, then the 
time-independence of the distribution of X is not enough because X involves, through its 



definition Eq. (16), a time- dependent transformation from the laboratory frame, in which 
the physical parameters are measured, so that they would still depend on time. Only when 
these parameters have become stationary can one say that the system is stationary. 

The SSFT was originally formulated for the average entropy production fluctuations on 
trajectory segments of length r along a single trajectory in the stationary state. Here we 
consider the statistics of the total work done on the system over time r, divided by /c_bT, 

Wr = P ^^'dtv^F(xi,x,*), (32) 
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for a sequence of initial times U of segments, all of length r, along a single stationary state 
trajectory = 1, 2, 3, . . .). To get the distribution of Wr of the segments along a trajectory, 
we use the following reasoning. According to the Eqs. (0) (for the force) and (|T6D (the 
definition of X), the expression in Eq. ( p2D is linear in Xj (just as in the transient case) and 
[as X( obeys the Langevin equation (|T^] we still have a Gaussian Green's function and a 
Gaussian stationary state, so that the distribution of Wr for each ti is again Gaussian: 



g 2Vt^ (r) 



PuiWr) = , (33) 

with the mean and the variance given by, respectively: 

M,Xt) = -k(3 r^^dt^l ■ (y* - x,*), (34) 
Jti 

VUr) = 2kf3 j'^^^^^dt', iyt[ v,*, ■ V,*, e-(*^-*'i)/-. (35) 

We assume that for sufficiently large ti, Mt- and Vt- will reach steady state values (see Sec. |T] 
for examples), and become independent of i. If in addition, the correlation between different 
segments {[ti,ti + r] and [tj,tj + r], say) decays sufficiently fast (when \ti — tj\ gets larger), 
then the distribution of Wr along a trajectory in the stationary state is given by: 

[Wr-Mg(-r)f 

PsiWr) = ^ . (36) 

Here the subscript S denotes that this distribution refers to the distribution of Wr over 
segments along the stationary state trajectory. The mean Ms is, from Eq. (0) and using 
Eq. (p!9D, given by: 

rt+T rt'2 

Ms{t)= \imk(3 / dt'2 / dt[e'^'''~''^^/^^v;, ■ v*, , (37) 
Jt Jo ' ' 

while the variance Vs is, from Eq. (|35D , given by 

Vs{t) = \im2k(3 dtU c/t'ie^(*2-*'i)/"' v*, • V*,. (38) 
Jt Jt ' ' 



Note that in the inner most integral in the expression for the mean in Eq. (pTj), the lower 
bound extends to time zero, whereas in the expression for the variance in Eq. (pSj), it extends 
to t. This is the origin of the fact that Vs and 2Ms are not identical [while Vr = 2Mt, 
Eq. (pQl)]. The deviation can be characterized by 

^ 2M,(r) - Vsir) 
'^^^ - 2Ms{t) ^^^^ 

Using this definition and Eq. (|36|) , one sees that 

PsiWr) \ Wr \ 



Ps{-Wr) ' 1 1 - e(r 
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This means that provided that 



£:(r) ^ as r ^ oo, (41) 

we have 



Vsir) 2Ms{t) as r ^ oo, (42) 

and the SSFT holds 

Pr.(W \ 

e^^ as r ^ oo. (43) 



Ps[-Wr) 

Of course, for any given v^, Eq. ( |4T| ) can be tested, but how general can we expect it to 
be satisfied? We write thereto Eq. ( pQ]) with Eqs. (|37|) and ( p8D as: 



e[T) = 



M.(t) 



hm fc/3(x:-y,*) f dt',e-''^l^^^l^, 

lir) ■ 

where Eq. ( |T9|) has been used. Here, the denominator is the total work done of the system 
in the stationary state in time r. If we are not in equilibrium, this is positive and grows 
with r. In the numerator, the exponential in the integral will make the integral bounded 
for large r, provided that does not grow exponentially in time with an exponent bigger 
than . Then e will become zero oc 1/r as r approaches infinity, and the SSFT in Eq. (^Sf) 
holds. 



D. Integrated Fluctuation Theorems 

In experiments such as done by Wang et a/. it is easier to check an integrated fluctua- 
tion theorem[|l^], because it is easier to obtain then good statistics for the required quantities. 
The Integrated Transient Fluctuation Theorem (ITFT) reads 



Pr(H^r < 0) 
Pr^r > 0) 



(45a) 



where the left hand side is the quotient of the probabilities to see a negative resp. a positive 
total work Wr after a time r: 



PTiWr<0)= / dWrPxiWr) (45b) 
J — oo 

PriWr > 0) = 1 - PriWr < 0) (45c) 
and the right hand side of Eq. (|45a|) is the average of exp{—Wr) over positive Wr, i.e., 

J,^dWrPT{Wr)e-'^^ 



j;;^dWrPTiWr 



(45d) 
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The ITFT of Eq. ( [45a| ) can be derived from the TFT in Eq. (|3l|) by first rewriting 
Pt{W^ < 0) = j\dWrPT{Wr) as 



/ dWr PT{-Wr)e^^ 
J — oo 

dWr PT{Wr)e~^\ 



(46) 



and then dividing by PriWr > 0). 

An Integrated Stationary State Fluctuation Theorem (ISSFT) can also be derived, but 
it is a little more subtle. Thereto, one has to consider whether 



PsjWr < 0) 
PsiWr > 0) 



(47a) 



holds, where 



Ps{Wr<0) ~ I dWrPs{Wr) 

J —CO 

Ps{Wr>^) = l-Ps{Wr<0), 



and 



f;;^dWrPs{Wr)e 



-Wr 



(47b) 
(47c) 

(47d) 



To start the derivation of the ISSFT of Eq. (|47a|) , the numerator of Eq. ( [47d|) is rewritten, 
using Eq. (^), as 



dWr Ps{Wr)e-^^ 
dWr Ps{-Wr) exp 



dWrPsiWr)exp 

[Wr-Ms{T)]'^ 



l-.(r) 

eir)Wr 
l-e{T) 



exp 



2Vs{t) 



e{r) 
l-e(r) 



(48) 



where Eq. ( P^ ) was used. We saw that the SSFT holds if £ — > for large r. Consider the 
exponent in Eq. (^8|). Writing out the square, this has a term linear in W^- of the form: 



Ms{t) 



Vs{r) l-e{r) 



Msjr) 
Vs(t) 



[l~2e{T)]Wr 



(49) 



As r ^ oo, we can neglect e compared to 1. Since this is the only place where e occurs, we 
can set e equal to zero on the right hand side of Eq. (|48|) , which then becomes P{Wr < 0). 
Dividing by P{Wr > 0) on both sides in Eq. (|4|) now yields the ISSFT in Eq. (|7|). 
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For the purpose of the investigation of the observabihty of the fluctuation theorems, 
discussed in Sec. |TTT|, we end this section by giving the exphcit forms of the left and right 
hand sides of the integrated fluctuation theorems Eq. ( [45a| ) and Eq. ( |47a| ). Defining 

^(^^-Wt^^ «.W.(e->;; (50a, 
^^M^i^^ «.M-(e->:. (50b) 

the TFT states that Lt = Rt, and the SSFT that Ls = Rs (the latter for large r only). 
Using Eqs. (|22|) , (|3^) and ([45 a]) , we get the following explicit expressions: 



Lt{t) = Rt{t) (51a) 



1 — erf 



Mt(t) 



Rt{t) = y (- (51b) 

1 + erf 



1 — erf 
1 + erf 



Ls{t) = ) f (51c) 

/ Msjr) 



1 — erf 
1 + erf 



Vsir)-Ms(r) 



Rsir) = e^--sir) \ V^)^. (51^) 

/ Ms{t) 
\yj2Vs{T)^ 



We can simplify the expressions for Lt and Rt using the relation between Mt and Vt in 
Eq. 

1-erf (iV^^) 

LAr) = ) ^ ( ■ (52) 

1 + erf [^l^/M^j 

In order to demonstrate the difference between Ls and Rs, we rewrite Eq. ( ^ld|) , using 
Eq. (p9D, in terms of e as 

[1-2s{t)]Ms{t) 



1 — erf 
1 + erf 



which shows that only for r — oo, Ls = Rs, i-e., that the work related ISSFT holds. 



E. Transient Fluctuations versus Stationary Fluctuations 

An interesting relation can be derived for the ratio of the probability of a negative total 
work and that of a positive one, for the transient case {Lt) and the stationary case {Ls). 
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Using Eqs. (|51al -c) and the asymptotic expansion of the error function, 



e 



erf (x) = 1 - ^ [x-^ + 0{x-^)] , (54) 

one obtains 




^ ' ^ ^.Q 2Vy(r) + 2Vs(r) 



Ls ^VsMt" ■ ^^^-^ 

Here the 0{x~'^) in Eq. could be neglected. This, because when r — oo, Mr and 
Ms [Eqs. (p^ and ( p7D respectively] will both grow linearly in time, Mt ~ Ms ~ 0{t) 
and similarly Vr ~ ~ 0{t) [by Eqs. (pQ]) and (|4^)], so that the arguments of the error 
functions in Eqs. ( p lb] ) and ( plc| ) become large (~ ^r) for large r. In fact, Ms and 
will grow with the same coefficient, as the rate of work done on the system will become 
stationary, but they will in general have a bounded difference, i.e., Ms — Mj- ~ as 
will Vs and Vr, i.e., — ~ One therefore has, using Eq. (|29|) and Eq. (^2]), an 

asymptotic behavior 

Mt{t) ^ wr + ai, Vt{t) ^ 2wt + 2ai 

Ms{t) ^ WT + a2; Vs{t) ^ 2wt + (56) 

for r — > oo, where the are independent of r. w is in fact the asymptotic rate at which work 
is supplied to the system, i.e., the consumed power. By expanding in terms of 1/r, we obtain 
for the exponent in Eq. (||) : -[M^{t)/2Vt{t)] + [MI{t)/[2Vs{t)] = -{ai + \a2 - \a^, i.e., 
it approaches a non-zero constant. The prefactor in Eq. (|55|) can similarly be shown to go 
to one as r oo, so that we have for large r. 



Lsir) 



Mg(r)-Mr(T) Vg(r)-Vr(r) 

e ^ 8 , (57) 



where Eq. (|S^) has been used to re-express this ratio in terms of the means and variances 
of the transient, resp. stationary state. Because the right hand side of Eq. (|57D is not equal 
to 1 in general, this shows that the ISSFT is not just the limit of the ITFT for large r. 



III. APPLICATIONS 

In this section, two kinds of motion are considered for the harmonic potential. In both 
cases parameters will be varied to see under what conditions an experiment would be able to 
demonstrate the integrated work fluctuation theorems (both transient and stationary) most 
convincingly. The motion that will be considered first, corresponds to the situation in the 
experiment of Wang et ai, i.e., it is a uniform linear motion. The other is a circular motion 
and might be implemented in a future experiment. Both approach a stationary state. 



A. Linear Motion of the Harmonic Potential 



The particular case considered here is a linearly moving harmonic potential, i.e., xj" = 
Voptt^ for t > 0. The quantities with which to test the work fiuctuation theorems, are given 
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in Eqs. (|51a| -d). The only unknowns are the means and variances of the transient and 



stationary distributions and these are given by Eqs. (PH), (]29|), (p^) and (pSl). If we insert 
Vj, which is a constant foptx here, into these equations, we obtain straightforwardly 

Mt{t) = w {t - rr[l - e~^/^^]} (58a) 

Vt{t) = 2Mt{t) (58b) 

Msir) = WT (58c) 

Vsir) = Vrir) = 2w {r - 7:,[1 - e"^/^"]} , (58d) 



where 



w = apv^pt, (58e) 

which can be interpreted according to Eq. ( ^8c|) as the rate of work delivered to the system 
per unit time. The equality between and Vr in (|58d| ) follows because the velocity of the 
center of the harmonic potential is constant, so that the integrands in Eqs. (p8|) and (|38D 
only depend on the difference ^2 ~ ^'i; and the shift over t in the definition of Vs is irrelevant. 
For this case, Eqs. ( |58a| -d) were already derived by Mazonka and Jarzynski in Ref. [jl2| 



By the theory presented in the Sec.|l|, the TFT holds for any motion of the harmonic 
potential, hence also in this case Lt = Rt- The SSFT holds if e [Eq. (pQ])] vanishes as 
r — *■ cxD, and this is so here, since 

T 

We now discuss the observability of the work related ITFT and the ISSFT for this model. 
There are only two relevant parameters for the fluctuation theorems here, the relaxation time 
Tr [Eq. (|13[)] and the rate of dimensionless work done w [Eq. ( |58e|) ]. To obtain realistic values 



for these parameters, orders of magnitude of various quantities can be taken from Ref. 
With the radius R of the order of three microns, and the viscosity of water rj of the order 
of 10~^kgm~^s~^, according to Eq. (|^), a is of the order of 5 x 10~^kg/s. With k of the 
order of 10~^kgs~^, Tr becomes of the order of 0.5s [Eq. (0)]. Furthermore, taking the 
temperature to be 300 K, gives (3 of the order of 2.4 x 10^°kgm^s~^, and with Vopt of the 
order of 1 nm/s, we find from Eq. ( |58eD that w is of the order of 12 s~^. 

For the case that w = 12 s^^ and r^. = 0.5 s, the expressions in Eqs. (|51a| -d) using 
Eqs. ( ^8a[ -d), are plotted together in Fig. |l|a. It is striking that even though we know that 
the ISSFT holds for sufficiently large r, this is not at all observed in the figure: the curves 
of Ls and Rs are completely different; in fact, the curve for Ls is indistinguishable from the 
T axis. Furthermore, both of these curves are different from Lt, which, given the result in 
Eq. (|57|) of Sec. |II E| , is less of a surprise. Clearly, the range of r for which the ISSFT is valid 



lies beyond the point where both Ls and Rs have relaxed to zero in Fig [l|a. This means 
that for the parameters typical of the Wang experiment, the ISSFT cannot be observed, in 
contrast to the ITFT, as the curve of Lt can be seen clearly, and Lt equals Rt, so that the 
ITFT could be observed. 

The reason that there is such a big difference in the signal for the transient and the 
stationary case, i.e., that the negative work fluctuations are more suppressed in Ls than 
in Lt, is the following. Since Mt{0) = [Eq. (P^l)] in the argument of the error functions 
in Eq. (0), Lt{t = 0) = 1. As the function Lt decays with increasing r, the question of 
whether it can be observed depends on whether it does not decay too quickly. On the other 
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FIG. 1: Integrated Fluctuation Theorems for the work for the linearly moving harmonic potential: 



Lt = Rt [cf. Eq. (51a)], Rs and Ls versus r (varying ranges) with (a) Tr = 0.5s and w = 12s' 



-1 



(b) Tr = 0.2 s and w = 4.8 s-\ (c) Tr = 0.08 s and w = 1.92 s-\ and (d) Tr = 0.032 s and 
w = 0.768 s~^ (in the last case the curves of Lt = Rt and Rs are indistinguishable). 



hand, the argument of the error function in Eq. ( |51c| ), Ms/ ^/2Vs does not have a limit of 
zero for r — > 0. In fact, using Eqs. ( |58c| ) and ( ^8d|) and expanding in r, one finds 



\/wTr/2, 

which can be large. So with Eq. ([51c|), one has for r ^ 0: 



Ms{r) 
lim — , 



(60) 



Ls{Q) 



l + eii{^JwTr/2)' 



(61) 



For the case plotted in Fig. |l|(a), ^JvotJ2 = ^/3, so that Ls{Q) = 7 x 10~^. No wonder we 
cannot see it in Fig. |l](a). Ls is exponentially suppressed for large values of WTr, which is 
the average work done in the stationary state during a relaxation time r^. As Wr = for 
r = for all trajectories, it follows from definition ( |50b|) that Rs{0) = 1- Comparing this 
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FIG. 2: Illustration of the relation between transient and stationary work fluctuations [Eq. (| 
for the linearly moving harmonic potential. Lt and are plotted logarithmically as a function of 
r, with Tr = 0.2 s and w = 4.8 s~^ [cf. Fig. |l]b]. The constant ratio in Eq. (57) becomes a constant 
distance between the curves in this logarithmic plot. 

value with the value of ^5(0) [Eq. (|6l|) ], we see that there is no chance to observe the SSFT 
if the work done during is too large. 

Thus in order to observe the SSFT we have to reduce the work done in the time Tr. One 



direct way to do this is to change the particle's radius. From Eqs. (H) and (0), we see that 
Tr oc R, whereas from Eq. ( |58e| ) , w oc i? as well, so WTr oc i?^. Another way would be to 
reduce the velocity of the harmonic potential, which does not affect r,., but changes w oc v^^^. 
Choosing to work with the radius as the control parameter (partly because smaller particles 
than used in the Wang experiment are commercially available), we plotted in Fig. |l|b-d, 
what happens when we make the particle 2.5 times smaller consecutively, thus reducing the 
work done per relaxation time each time by a factor 6.25. In first instance [Figs. |I|b and |l]c] 
the curve of Ls gets closer to the curve Rt, and Ls starts to get visible. But only in the last 
graph. Fig. |l]d, where the particle's diameter is about 15 times smaller than in Fig. |l|a (which 
would mean R ^ 200 nm in the Wang experiment), do we clearly see that Ls approaches 
Rs- In this case, WTr ~ 0.02, confirming that the work done in needs to be small to see 
a convincing signal of the ISSFT. 



Finally, we look at Eq. (57), which says that in the long r limit, the ratio of Lt and Ls 



becomes a constant. Using Eqs. (|58a| -d) in this case, Eq. (|57D reads 

Taking = 0.2 s and w = 4.8 s~^, we plotted L^ and Ls logarithmically to illustrate 
Eq. (0). Eq. (|62D shows once more that negative fluctuations of the work are more sup- 
pressed in the stationary state than in the transient period. 
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B. Circular Motion of the Harmonic Potential 



In the case of a circular motion of the harmonic potential, we write 

X* = r{sin(nt)x + [1 - cos((]t)]y}, (63) 
for t >0. We use the same procedure as in the previous case, i.e., we determine 

V* = rQ{cos{Qt)ii + sin(^]^)y}, (64) 
for t > 0, insert this into Eqs. (0), (0) and (|55|), and obtain straightforwardly 



, , , 2^]r,. sin(fir)e-^/^' 

Mt{t) = w{t - Tr- ^ ' 



where now 



1 + nH^ 

[l-^^V][l-cos(r]r)e-/-1 1 

TTT^^^f / ^^^^^ 

Vt{t) = 2Mt{t) (65b) 

Msir) = WT (65c) 

Vs{t) = Vrir) = 2MT(r), (65d) 



1 + fi2r2 

The expression for e follows using its definition Eq. (^ 



.(.) ^ ^^^^|21^..sm(r^r)e-^^ 

+ [1 - n^T^][l - cos(r]r)e-"/"'-]| (66) 

which again vanishes like 1/r when r oo, so that the ISSFT holds, i.e., Ls = Rs for large 
r. 

Note that the Eqs. (|65a| -e) and (|^) reproduce the Eqs. (|58a| -e) and (|59|) in the limit for 



r2 ^ 0, keeping Vopt = rQ constant. But these equations are not just an extension of the 
linear case. Under the resonance condition fir^ = 1, e in Eq. (|66|) becomes 

e{T) = iL^ilMe-r/rr^ (67) 
r 

i.e., it decays exponentially, rather than oc 1/r (which it does for all other choices of Q). In 
addition, in the resonance case e is zero at times mrTj. for n = 1, 2, 3, . . ., and exponentially 
small (~ e""''") in between; consequently, at those times, Ls and Rs are equal [Eqs. ( ^lc|) 
and (|53|)], whereas they are exponentially close in between. This means that the ISSFT 
holds and can be visible at much shorter time scales than in the linear case. 

In Fig. ^ Lt = Rt, Ls and Rs from Eqs.( |51a[ -d) are shown for two cases, which for 
= become identical to the cases (6) and (c) in Fig. |l|. For these cases we show what 
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FIG. 3: Integrated Fluctuation Theorems for the work for the circular motion Lt = Rt [cf. 
Eq. ( plap ], Ls and Rs versus r with (b') = 0.2 s, w = 4.8 s~^ (cf. Fig |^b) at resonance 
= l/r^, (b") same and w but Q = S/r^, (c') = 0.08s, w = 1.92 s~^ (cf. Fig. ||c, Rs is 
indistinguishable from Lt = Rt) at resonance Q = l/Tr-, and (c") same and w but 17 = 5/rj.. 
Note that in (b') and (c'), beyond a time TTTr, all curves become indistinguishable. 

happens at their respective resonance points Vt = l/r^ and what happens when VL is larger, 
5/r,.. In varying the frequency, we keep w fixed, which physically means we would have to 
adjust r according to Eq. ( |65e|) . 

From Fig. ^d', we see that at the resonant point VL = l/r^, while Ls could not be seen 
in the f» case (cf. Fig. |l|b), it can now be seen, but is still too small compared to Rs 
to check the SSFT. If at = the value of Ls could be seen (Fig. |l]c), then going to the 
resonance, improves the situation: in Fig. ||c', Ls and -Rs approach each other after a far 
shorter time than in Fig. |I]c, and become indistinguishable after r = ttt^, because of the 
exponentially small difference between them as discussed above. We remark that when the 
agreement is already good in the f2 = case (like in Fig. |l]d), going to the resonance changes 
little. Furthermore, we see from Fig. ^ that going beyond the resonance {Q = S/r^) will 
cause the curves to deviate from each other again. We also note that above resonance, the 
work fluctuations in the stationary state are larger than those in the transient case, which 
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is the opposite as for the hnear motion. In fact, the form Eq. ([57| ) takes here is [using 
Eqs. (1^), (|65^) and (|65H|) ] 



for r — > oo. The exponent changes sign going from Qtj. < 1 to Qrr > 1. What is hap- 
pening physically is that the system is driven so fast that it cannot relax within one cycle. 
This increases the fluctuations Wr- On the basis of the foregoing, we see that in order to 
demonstrate the ISSFT, it helps to go to a resonant circular motion. 



IV. DISCUSSION 

1. Inspired by the experiment of Wang et al.^ which showed that the work related 
ITFT holds when a small latex bead is dragged linearly through a fluid by means of a laser- 
induced harmonic force, a model of a Brownian particle in a harmonic potential with an 
arbitrarily moving minimum was used to study both the work related ITFT and ISSFT, 
under more general conditions. From the Langevin equation that describes the motion of 
the Brownian particle in this simple model, everything can be explicitly calculated. As 
expected, the work related TFT holds for all time, as does its integrated variant (ITFT). 
The work related SSFT and its integrated version (ISSFT) hold for sufficiently large times 
provided a stationary state exists. 

We also found a new relation between work related ISSFT and ITFT. If one looks at the 
ratios for the probability to find in a time r a negative vs. a positive work done on the 
system, in the transient state {Lt) and in the stationary state {Ls), then Lt/Ls approaches 
a constant (which is not 1) as r — > oo, given by Eq. dSTf ). 

2. We have not found many choices for the motion of the harmonic potential (xj') for 
which a stationary state exists [in the sense that the limits of Eqs. (^) and (^) exist]. 
There is the linear motion corresponding to the Wang et al. experiment, which has been 
worked out in Sec. |III A|, and there is the possibility of a circular motion treated in Sec. |IIIB 



as well as a spiral motion, which is a trivial superposition of the previous two. However, any 
motion which in the course of time approaches one of these cases, will also reach a stationary 
state, corresponding to that case. In contrast to the simple motion of the harmonic potential 
considered in Sec. |T|, allowing for an arbitrary motion of the harmonic potential, may give 
rise to (arbitrarily) different fluctuations in the transient and the stationary cases. 

3. We note that on the basis of our explicit calculations, we are able to explore under 
what conditions the work related ITFT and ISSFT might be observable, which is relevant 
for the devise of future experiments. 

For the ITFT, which holds for all time, observability is purely a matter of how fast the 
quantity Lt decays: if it decays too fast, the ITFT will not be measurable. We showed 
that if one inserts values taken from the Wang experiment [^ into the explicit expression of 
Lti the ITFT shows a clear signal which decays on the order of a second, consistent with 
the fact that the ITFT could be observed in that experiment. The relaxation time, i.e. t^., 
however is off; the relaxation time found in the experiment is of the order of 1 to 2 seconds, 
whereas the value of the harmonic force constant and the application of Stokes' Law give a 
relaxation time of 0.5 seconds. This could be due to boundary effects, local heating due to 
the laser, deviations from the harmonic nature of the laser-induced force. 
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The fact that the ITFT can be observed does not imply that the ISSFT can be observed. 
In fact, problems with observability arise when we insert values taken from the Wang ex- 
periment into the quantities for the ISSFT: the signal of Ls is then too small because 
the work done on the system in a relaxation time is too large. There are a few ways to 
improve this situation, i.e., to make the ISSFT observable in an experiment: first, one can 
take a smaller particle, or second, make the velocity with which it moves through the fluid 
smaller. Both methods reduce the work done in a relaxation time r^. If the diameter of the 
particle is reduced to about 16 times smaller than the original one (which means about 400 
nm across), we see that the ISSFT can indeed be observed. 

The circular motion that we investigated offers a third possibility to improve the observ- 
ability of the ISSFT: Under resonance conditions [VLTr = 1), the deviations from the ISSFT 
become exponentially small after a time Trr^. 

4. We end by giving some issues which are open for future investigation. Some possi- 
ble extensions of the theory could be the following. The theory developed here is for the 
overdamped case only. One might wonder if there is ever any practical need to consider the 
situation where the damping is not so large. That would mean the theory would start with a 
Langevin equation of motion for x and v, or a Kramers equation. We have carried out such 
calculations for the case of a linear motion of the harmonic potential and found the same 
results as reported here. In the case of the circular motion, this would also allow a more pre- 
cise discussion of the role of the centrifugal force (due to Q), which a rough estimate limits 
to Q -C a/m (this is of the order of lO^Hz for parameters taken from the Wang et al. ex- 
periment). One could also consider the case of anharmonic rather than harmonic potentials. 
On a practical level, is seems plausible that the theory could be applied to other systems, 
such as linear electrical circuits and could be particularly relevant to nano-technology. 

On a more fundamental level, one could ask what the precise relation is of the work 
fluctuation theorem (for PW^°^) discussed here and the usual entropy production theorems 
(for /5iy^^) for dynamical and stochastic systems[0, ^, |§, |§, 0]. While the work fluctuation 
theorems hold for all classes of systems considered so far, this appears not to be the case for 
the entropy production theorems. In a future publication, we intend to discuss this question 
in detail, since the theory needed for this deviates too much from the present one to be 
included it in this paper. We can state, however, that for the models considered here, such 
a theory indicates that while an SSFT for the entropy production, i.e., for PW^"^ , appears 
to hold for long times as usual, the TFT for the entropy production seems to hold for long 
times only as well[^, and not as an identity for all times, as would be expected |[I9| . 

This work has been supported by the U.S. Department of Energy, under grant No. DE- 
FG-02-88-ER13847. 
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